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Abstract. I review recent selected developments in the theory and modeling of ultrarelativistic 
heavy-ion collisions. I explain why relativistic viscous hydrodynamics is now used to model 
the expansion of the matter formed in these collisions. I give examples of first quantitative 
predictions, and I discuss remaining open questions associated with the description of the freeze- 
out process. I argue that while the expansion process is now well understood, our knowledge 
of initial conditions is still poor. Recent analyses of two-particle correlations have revealed 
fine structures known as ridge and shoulder, which extend over a long range in rapidity. These 
correlations are thought to originate from initial state fluctuations, whose modeling is still crude. 
I discuss triangular flow, a simple mechanism recently put forward, through which fluctuations 
generate the observed correlation pattern. 



1. Introduction 

The Relativistic Heavy Ion Collider (RHIC) at Brookhaven has been running for ten years, and 
the Large Hadron Collider (LHC) at CERN is expected to collide nuclei before the end of 2010. 
In the three years since the last INPC conference, we have witnessed significant progress in our 
understanding of RHIC experiments, both at the quantitative and qualitative level. A recent 
breakthrough is that quantitative predictions from viscous relativistic hydrodynamics are now 
available [1, 2], which allow for a better description of the bulk of observables. Section 2 focuses 
on this topic, although there have been important developments in ideal hydrodynamics as well 
[3, 4, 5]. On the experimental side, new analyses involve rare particles, such as heavy flavours, 
and detailed investigations of the fine structure of correlations. I concentrate on the latter topic 
in section 3, and I discuss a new mechanism which sheds new light on the "ridge" and "shoulder" 
structures observed in correlations, which have puzzled theorists for a few years. 

2. Relativistic viscous hydrodynamics 

Nuclei colliding at RHIC are thin pancakes: the Lorentz contraction factor is 100 at the top 
energy. The clear separation between the small longitudinal size and the much larger transverse 
size R is a key feature of heavy-ion phenomenology. In a time short compared with R (in natural 
units c = 1), matter is produced: a system of strongly interacting quarks and gluons, which 
rapidly expands into the vacuum. The underlying microscopic description is a quantum field 
theory, quantum chromodynamics, in the nonperturbative regime, which does not necessarily 
admit a kinetic description in terms of quasiparticles. However, if the system is large enough, 
hydrodynamics gives a reliable description of its real-time macroscopic evolution, once initial 
conditions are specified, at a time t <C R. 



For a nonrelativistic fluid, the usual form of the Navier-Stokes equation is 



p^— + v- Vvj = -VP + 7]V 2 v. (1) 

where p is the mass density, v the fluid velocity, P the local pressure, and r] the shear viscosity. 
Let us analyze the order of magnitude of the various terms. Our analysis will differ from 
that usually found in textbooks. Usually, the Navier-Stokes equation is used in the context of 
incompressible flows. Incompressible does not refer to an intrinsic property of matter. It is a 
synonym of slow: it characterizes a flow where v is much smaller than microscopic (thermal) 
velocities v t h- In this case (which applies for example to the lower atmosphere), it can be shown 
that the density is almost constant throughout the fluid, hence the term "incompressible". In 
a heavy-ion collision, however, the fluid is freely expanding into the vacuum, which is a clear 
instance of a compressible flow. Transverse expansion is driven by the thermal motion, so that 
v ~ i>th- Gradients and time derivatives all scale like 1/R. In Eq. (1), all terms scale like 1/R, 
except the viscous term which scales like 1 / R? . Its magnitude relative to the acceleration term 
is X/R, where A is a length involving thermodynamic quantities. Inspection of equation (1) gives 
A = rj/pvth- In systems which can be described by kinetic theory, A is the particle mean free 
path. Now, fluid dynamics is a macroscopic description which is valid only if A <C R- Therefore 
the viscous term must be a small correction, or fluid dynamics itself breaks down. 

The dynamics of compressible fluids can be formulated as a systematic gradient expansion, 
that is, an expansion in the Knudsen number K = X/R [6]. Ideal hydrodynamics is the leading 
order, and Navier-Stokes theory is the first correction of order K. The second-order correction, of 
order K 2 , is usually called Israel-Stewart theory in the relativistic case. The expansion converges 
if K is small: hydrodynamics applies to heavy-ion collisions only if the nucleus is large enough, 
and if the viscosity rj is small enough. 

Until 2007, the state of the art for RHIC phenomenology was ideal hydrodynamics [7], 
motivated by the hope that the viscosity of hot QCD is very small [8, 9]. It was realized early 
on that viscous corrections were likely to be important, both on theoretical grounds [10] and 
from a critical analysis of RHIC data [11]. However, going beyond ideal hydrodynamics proved 
a difficult task. Relativistic Navier-Stokes equations are known to break causality [12] because 
small-wavelength modes propagate faster than light. This is not a conceptual problem, since 
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Figure 1. (courtesy I. 
Bouras [13]) Compari- 
son between second-order 
viscous hydrodynamics 
(vSHASTA) and kinetic 
theory (BAMPS), for two 
values of the viscosity to 
entropy ratio rj/s. The 
dashed line is the initial 
condition at t = 0, and the 
curves display the pressure 
profile at t = 3.2 fm/c. 



hydrodynamics is only expected to describe long- wavelength modes [14], but it is a practical 
difficulty because of the associated numerical instabilities. As a consequence, one must include 
second-order terms. There were early attempts to apply second-order viscous hydrodynamics to 
heavy-ion collisions [15], but full calculations, including comparison with RHIC data, became 
available only recently [16, 17, 18, 19, 20]. There is some arbitrariness in second-order terms, but 
the goal is to estimate quantitatively the first-order correction, in regimes where second-order 
terms are negligible. 

If second-order terms could be determined from first principles, they would extend the validity 
of the hydrodynamic description up to larger values of the Knudsen number. Although this is 
not yet the case for QCD, a feasibility study was recently performed [13] for a dilute system 
described by kinetic theory, where second-order terms were determined unambiguously. Figure 1 
displays the comparison between hydrodynamics (labeled vSHASTA) and kinetic theory (labeled 
BAMPS). At time t = 0, the fluid is at rest, and its pressure is uniform with a discontinuity at 
z = (shock tube problem). At t > 0, the fluid flows toward lower pressures, i.e., to the right. 
In ideal hydrodynamics, the flow profile consists of three distinct parts: a continuous rarefaction 
wave propagating to the left, a compression shock (pressure discontinuity) propagating to the 
right, and a region of uniform flow in between. This pattern is clearly seen with the smaller 
value of the viscosity rj/s = 0.01. If the viscosity is increased by a factor of 10, the flow pattern 
is smoothed. In particular, the width of the shock is proportional to the viscosity. In both 
cases, there is almost perfect agreement between hydrodynamics and kinetic theory. Agreement 
is found to hold within 10% up to K = 1/2. Both the time scale in figure 1 and the larger value 
r]/s = 0.1 are of the order of expected values at RHIC, so that the difference between rj/s = 0.01 
and rj/s = 0.1 results gives an idea of the magnitude of viscous effects at RHIC. 




Figure 2. Average multiplicity of pions, 
kaons, protons, and antiprotons per partic- 
ipant, versus number of participants (from 
[1]). The number of participant nucleons is 
used to estimate the centrality of the collision. 
It is maximum for central, head-on collisions. 



Figure 3. Average transverse momentum 
of particles versus number of participants 
(from [1]). As in figure 2, experimental data 
from PHENIX [21] are compared with viscous 
hydrodynamics calculations. The label CGC 
refers to a specific choice of the initial density 
profile in hydrodynamics. 



I now discuss hydrodynamic calculations for collisions at RHIC. The initial energy profile in 
the transverse plane cannot be determined from first principles because interactions are largely 
nonperturbative. This is a major source of uncertainty. Most calculations use a simple ansatz 
which reproduces the observed centrality dependence of the pion multiplicity (figure 2). The 



multiplicity scales essentially like the number of nucleons participating in the collision. Particle 
ratios (proton to pion and kaon to pion) are also essentially constant. In hydrodynamics, they 
are determined by the freeze-out temperature at which the fluid transforms into hadrons. The 
average transverse momentum (pt) of identified particles (figure 3) also shows little centrality 
dependence. This observable reflects the kinetic energy per particle. In hydrodynamics, the 
kinetic energy is proportional to the effective temperature of the system, defined as the average 
temperature at a time of order R when transverse expansion fully develops. This temperature 
is determined (through the equation of state) by the density of particles per unit volume, 
which is itself closely related to the multiplicity per participant plotted in figure 2. Since this 
density depends little on centrality, hydrodynamics naturally predicts that (pt) depends little on 
centrality [22]. While the absolute scale of {pt) depends little on viscosity, it depends crucially 
on the equation of state used in the hydrodynamic calculation, which is taken from lattice QCD 
calculations. The equation of state determines the temperature as a function of the density, and 
the fact that hydrodynamics reproduces data is a non trivial test of the equation of state from 
lattice QCD. 



Figure 4. Elliptic flow scaled by 
the initial eccentricity e and by the 
hydrodynamic limit h (see text) as 
a function of the density of parti- 
cles per unit area. Results from 
ideal and viscous hydrodynamics [1] 
are plotted together with data from 
RHIC (from [23]). Open symbols 
and full symbols correspond to two 
different choices of initial condi- 
tions (called CGC and Glauber). 
Dashed and solid lines are fits us- 
ing equation (3). 

I now discuss an observable which is more sensitive to viscosity, namely, elliptic flow. The 
azimuthal distribution of emitted particles is elliptic, due to the elliptic shape of the interaction 
region in the transverse plane (see figure 9). The fluid velocity is proportional to the pressure 
gradient, which is larger along the smaller dimension. Particles tend to be emitted parallel 
to the fluid velocity, which results in more particles emitted along the minor axis [24]. The 
corresponding observable is 

v 2 = (cos(20-2^ 2 )), (2) 

where ^2 is the minor axis of the participant ellipse in figure 9, <p is the azimuthal angle of 
an outgoing particle, and angular brackets denote an average over particles and events. i>2 is 
smallest for central collisions, where the overlap area is circular except for fluctuations, and 
increases up to 7% for more peripheral collisions. Since elliptic flow is generated by pressure, it 
is a clean signature of collective flow. 

Figure 4 displays the centrality dependence of elliptic flow, plotted in such a way that results 
can be easily understood by dimensional arguments. The vertical axis displays t>2 scaled by the 
initial eccentricity, e, of the participant ellipse in figure 9. The horizontal axis is the multiplicity 
of produced particles scaled by the area of the ellipse. This density per unit area is often used as 
an alternative measure of the centrality. Since the effective density of particles per unit volume 




depends little on centrality, {1/ S){dN / dy) is proportional to the system size R: the largest value 
corresponds to central collisions, as in figures 2 and 3. In hydrodynamics, v^/e depends little on 
the details of initial conditions. It is essentially determined by the size R and by thermodynamic 
quantities, such as the speed of sound and the viscosity. Thermodynamic quantities depend on 
the temperature, but we have seen that the effective temperature depends little on centrality, 
so that V2/E depends on centrality only through the size R. In ideal hydrodynamics (top line 
in figure 4), v^js is independent of the system size. This is a natural consequence of the 
scale invariance of equation (1) with 77 = 0: velocity and pressure are unchanged under the 
transformation [t,x) — > (at, ax), where a is arbitrary. With a nonzero viscosity, thermalization 
is incomplete. Therefore, there is less pressure and less flow, and V2/E is smaller. Using the 
dimensional analysis explained at the beginning of this section, the viscous correction is expected 
to scale with r\ and R like i]/R oc r/((l/ S)(dN/dy))~ 1 . A look at the viscous hydro results for 
r]/s = 0.08 and rj/s = 0.16 in figure 4 shows that the scaling is at least qualitatively correct: 
deviations from ideal hydrodynamics are proportional to the viscosity, and inversely proportional 
to (l/S)(dN/dy). In order to check the scaling quantitatively, each set of results (given value of 
rj/s and choice of initial conditions) is fit using the two parameter formula 



The parameter h is the ideal hydro limit. The vertical axis in figure 4 is scaled by h in order 
to single out the viscous correction, quantified by the parameter B. This parameter is found to 
be proportional to rj/s, as expected. Unlike (pt), elliptic flow is very sensitive to the viscosity. 
rj/s = 0.08 is close to an absolute lower bound [9], yet it produces sizable effects, because an 
atomic nucleus is not quite a macroscopic object. Note that both sets of initial conditions (CGC 
and Glauber) give almost the same v-^/e, although V2 and e differ. 

The next step is to put experimental data on this plot. V2 and dN/dy are measured, but e 
and S must be calculated within a model, which entails uncertainties [25]. Two different models 
of initial conditions are used, and the results are then fit using equation (3). Both models take 
into account eccentricity fluctuations [26] due to the finite number of nucleons (see figure 9 and 
section 3), which are important for central collisions. Equation (3) fits experimental data quite 
well [27]. By comparing theory with data, we can extract the viscosity of QCD from this plot. 
Depending on which set of initial conditions is used, the value of rj/s is 0.16 or significantly larger. 
Our poor knowledge of the initial density profile hinders our ability to extract thermodynamic 
quantities from the data. With CGC initial conditions, the viscous correction is roughly 20% for 
central collisions. It is small enough for hydrodynamics to be a valid approach; but it is sizable, 
and viscosity must be taken into account in order to achieve quantitative agreement with data, 
in particular for elliptic flow. 

I conclude this section by mentioning some of the pending issues in hydrodynamics. While 
viscous hydrodynamics provides us with a detailed picture of the expansion, we still miss the 
beginning and the end of the story: the uncertainty in the initial density profile has already 
been pointed out. There are also uncertainties in the final state, when the fluid transforms into 
particles. This freeze-out procedure is well defined in ideal hydrodynamics: thermal equilibrium 
is achieved locally, so that momentum distributions are thermal at any point in the fluid rest 
frame. Viscosity leads to a small deviation from thermal equilibrium, and the momentum 
distribution in the fluid rest frame is 
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where I have neglected quantum statistics for simplicity. The dependence of 5 /viscous (p) on 
the direction of p is constrained by the continuity of the pressure tensor at freeze-out, but its 
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Figure 5. Elliptic flow versus transverse 
momentum in mid-central Au-Au collisions at 
RHIC. Calculations are shown for pions [28] 
using ideal hydrodynamics (77/s = 0.0001) 
and viscous hydrodynamics (rj/s = 0.16) with 
several choices of ^/viscous are shown together 
with recent data for charged hadrons from the 
PHENIX collaboration [29]. 
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Figure 6. The ratio of quadrangular flow, 
V4, to elliptic flow squared, versus transverse 
momentum. Theory as in figure 5. PHENIX 
data have been scaled down by a constant 
factor to take into account fluctuations of the 
initial geometry (eccentricity fluctuations) 
which increase U4/(t>2) 2 [30]. 



dependence on p = \p\ is not. Most calculations so far assume <5 /viscous oc P 2 [10], but it has 
been recently emphasized [31] that other choices are possible. The correct choice depends on 
the detailed structure of hadronic cross sections at freeze-out. Figure 5 displays a comparison 
between theory and data for elliptic flow. While ideal hydrodynamics predicts too large a 1)2, 
as expected from the previous discussion, viscous hydrodynamics matches the data. Results up 
to pt ~ 2 GeV are in fact fairly insensitive to the choice of <5/ v iscous> which is good news. For 
momenta larger than 2 GeV, results depend on <5 /Viscous, but the viscous correction becomes 
large and hydrodynamics breaks down anyway. Now, look at figure 6, which displays results for 
the fourth harmonic of the azimuthal distribution, is defined by 

v A = (cos(4<£ - 4* 2 )}- (5) 

Once scaled by (t>2) 2 , data are essentially flat as a function of transverse momentum. While ideal 
hydrodynamics naturally reproduces this flat behaviour, taking into account viscosity makes the 
agreement with data worse, even for moderate values of pt where hydrodynamics is expected to 
be reliable. Among the three choices of <5 /viscous) the standard quadratic ansatz, which is used 
essentially by all viscous hydro calculations so far, gives the worst results. 

Another shortcoming of present viscous hydrodynamics calculations is that they do not 
incorporate chemical freeze-out, i.e., the observation that ratios of particle abundances are 
constant below some temperature (typically 165 MeV). In ideal hydrodynamics, entropy is 
conserved, so that the entropy per particle is uniform throughout after chemical freeze-out. 
This makes chemical freeze-out easy to implement [32]. In viscous hydrodynamics, diffusion 
processes must be taken into account, resulting in complicated equations [33]. An alternative 
way of dealing with the hadronic phase is to couple hydrodynamics to kinetic theory below some 
temperature. This was done in the early days of RHIC for ideal hydrodynamics [34, 35, 36, 37, 38] 
and is being generalized to viscous hydrodynamics. This approach also has problems, because 
the coupling is done at a high density where the validity of kinetic theory is questionable. 
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Figure 7. Measured two-particle correlation 
in proton-proton collisions at centre-of-mass 
energy 200 GeV. [39]. 



Figure 8. Measured two-particle correlation 
in central Au-Au collisions, same energy per 
nucleon as in figure 7. 



3. Correlations and fluctuations 

The most striking differences between proton-proton and nucleus-nucleus collisions are probably 
found in correlations. I first define a measure of the correlation. In a given event, let Ni and 
N2 denote the numbers of particles in two different phase-space bins, say, with azimuthal angles 
4>\ and 4>2, up to some 5(f). The correlation is defined by 



C 



(N!N 2 ) 
(Ni)(N 2 ) 



(6) 



where angular brackets denote an average over a larger number of events within a centrality 
class. Correlations typically vary like the inverse of the system size: if a nucleus-nucleus 
was a superposition of n independent nucleon-nucleon collisions, C would be proportional to 
1/n. In order to remove this trivial size dependence, the PHOBOS collaboration multiplies the 
correlation by the average number of particles [40]. The correlation is measured as a function 
of the relative azimuthal angle A(f> = 4>i — (f>2, and of the relative pseudorapidity Arj. The 
pseudorapidity is a function of the longitudinal velocity, parallel to the beam. Results are 
displayed for proton-proton collisions in figure 7, and for Au-Au collisions in figure 8. The 
difference is striking. In proton-proton collisions, the correlation is peaked around Ar] = 0: 
correlated particles have similar longitudinal velocities. This correlation is essentially flat in 
azimuth, except for a small peak near A(f> = and a wider peak near A(f) = tt. In Au-Au 
collisions, on the contrary, the correlation has a marked structure in A(f), which extends to large 
Arj. 

Simple causality arguments [41] show that correlations at large Arj can only be created at 
early times. It is in fact likely that their source is contained in the wavefunction of the incoming 
projectiles. The A(f> dependence of the correlation function at large Ar] has two maxima at 
A(f> = and A(f> = tt, which can be explained by a term proportional to cos2A0. This term is 
due to elliptic flow. As explained above, the origin of elliptic flow lies in the elliptic shape of the 
overlap area. Assuming symmetry with respect to the direction of ifj2, one can rewrite equation 



(2) in complex notation, for a single event, 



{e 2i *)= V2 e 2l *\ (7) 

Assuming that the only correlation is due to flow, 

(cos2A</>> = (eMfo-fc)) = (e 2 ^i)( e - 2i ^) = (v 2 f . (8) 

This well-known correlation is usually subtracted from the observed correlation. 

A closer scrutiny of figure 8 reveals that the correlation at large Ar] is not exactly a cos 2 Acj) 
modulation. While the maxima at Acj> = and Acj> = it have approximately the same height, 
the distribution is narrower around the first maximum at Acf) = 0. The narrower peak around 
Acj) = has been dubbed the ridge and has puzzled theorists for several years [41, 42, 43, 44, 45]. 
One easily checks that this pattern can be reproduced by assuming, in addition to the cos 2A<p 
term, a term acos3A0 — 6 cos Acj) (see figure 11). The maxima at Acj) = and Acj> = tt have 
the same height if a = b. Transverse momentum conservation naturally generates a — cos Acp 
term [46], but there was no known mechanism to produce a cos3A0 correlation. 




Figure 9. Interaction region projected on 
the plane perpendicular to the collision axis. 
Circles indicate the positions of nucleons 
within nuclei just before the collision. Purple 
nucleons are the participants, which undergo 
at least one collision. The positions of 
participants define an ellipse marked as a red 
line. ^2 is the azimuthal angle of the minor 
axis, along which elliptic flow develops. 



Figure 10. (courtesy B. Alver) Another 
event, which has been chosen because the 
distribution of participant nucleons draws a 
triangle in the transverse plane. By analogy 
with figure 9, one denotes by ^3 the flat 
axis of the triangle. 



The breakthrough came a few months ago, when Alver and Roland [47] proposed such a 
mechanism. This mechanism involves fluctuations, which I have already mentioned, and now 
explain in more detail. Due to Lorentz contraction, the time scale of a heavy-ion collision at 
RHIC is so short that the positions of nucleons within nuclei are frozen during the collision. The 
collision acts as a quantum measurement process: the positions of nucleons are random numbers, 
whose probability is determined by the wave function of the nucleus. Now, each nucleon-nucleon 
collision produces several particles whose rapidity span a large interval. One therefore expects 
that at all rapidities, the initial density profile has "hot spots" at points in the transverse plane 
where there are more nucleons. These fluctuations are playing an increasingly important role in 
heavy- ion phenomenology. The new idea of Alver and Roland is that in some events, fluctuations 
may create a triangular shape (see figure 10). In the same way as elliptic flow develops along the 



minor axis of the ellipse, triangular flow develops along the flat side of the triangle. By analogy 
with equations (7) and (8), we write 

(e 3 ^) = V3 e 3i * 3 (cos3Ac/>> = (v 3 f. (9) 

In the same way as Vi is proportional to the initial eccentricity, turns out to be proportional to 
the triangularity, £3 [47], a dimensionless number which characterizes how triangular the initial 
distribution is. This triangularity is produced by statistical fluctuations. For central collisions 
shown in figure 8, however, the eccentricity also results from statistical fluctuations [26], and is 
of the same order of magnitude. One therefore expects a similar magnitude for both cos(2A</>) 
and cos(3A0) components. The triangular term is smaller because viscous damping is larger for 
V3 than for V2 [48]. However, ^3 increases more rapidly with pt than V2, and naturally produces 
a shoulder structure of the correlation function (see figure 11) with a high-p^ trigger, similar 
to the one observed experimentally [49, 50]. Note that triangular flow appears naturally in 
hydrodynamic simulations with fluctuating initial geometries [45, 51]. 
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Figure 11. Pair azimuthal correlation 
R(A4>) oc v|cos2A0+?;|cos3A0-6cosA0 
from elliptic and triangular flow. I have 
chosen b = \v\ for both curves. Shoulders 
appear around Acj) = ir if v% > \v2- 



4. Conclusions and perspectives 

There is now a wide consensus that heavy-ion collisions at ultrarelativistic energies produce a 
fluid with a small viscosity over entropy ratio. Relativistic viscous hydrodynamics describes 
the expansion of this fluid, therefore it should explain the bulk of the observables: momentum 
spectra, anisotropies of identified particles in the soft sector (below p t ~ 2 — 3 GeV) , as well as 
Bose-Einstein correlations (not covered in this talk) . The recent developments on triangular flow 
have shown that even the detailed structure of correlations is naturally explained by flow. At 
LHC, the higher density will result in an increased lifetime for the fluid, and heavy-ion collisions 
should be dominated by collective flow. On the theoretical side, viscous hydrodynamics is still 
in its infancy, and much work remains to be done in order to properly incorporate chemical [33] 
and thermal [28] freeze-out. 

Our understanding of nucleus-nucleus collisions is hindered by our poor knowledge of the 
initial state. In particular, quantum fluctuations, resulting from the fact that colliding nuclei 
are quantum objects, are important for the little bang. In fact, the role of fluctuations is in 
several respects the same as in big-bang physics [52, 53]. Present models of fluctuations are 
rather crude. There has been huge progress in our understanding of the wavefunction of a 
nucleus at high-energy, and this should be exploited fully to estimate quantities relevant to 
heavy- ion phenomenology, such as the initial eccentricity [54] and fluctuations. 
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